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Abstract 

The solutions of a renormalized BCS equation are studied in three space dimensions 
in s, p and d waves for finite-range separable potentials in the weak to medium coupling 
region. In the weak-coupling limit, the present BCS model yields a small coherence length 
£ and a large critical temperature, T c , appropriate for some high-T c materials. The BCS 
gap, T c , £ and specific heat C S (T C ) as a function of zero-temperature condensation energy 
are found to exhibit potential-independent universal scalings. The entropy, specific heat, 
spin susceptibility and penetration depth as a function of temperature exhibit universal 
scaling below T c in p and d waves. 

PACS Numbers: 74.20.Fg, 74.72.-h 

1 Introduction 

At low temperature, a collection of weakly interacting electrons spontaneously form large over- 
lapping Cooper pairs |1| according to the microscopic Bardeen-Cooper-Schreiffer (BCS) theory 
of superconductivity P, |3[]. There has been renewed interest in this problem with the discovery 
of enhanced superconductivity in alkali- metal-fulleride compounds 0] and cuprates ||. The 
fulleride compounds, with critical temperature T c up to ~ 30 — 40K, exhibit superconductivity 
in three dimensions and have a relatively small coherence length £ : ^kp ~ 10 — 100, with kp 
the Fermi momentum. At zero temperature £ is essentially the pair radius. In the application 
of the BCS theory to high-T c materials, the serious challange is to consistently produce a large 
T c and a small £ in the weak-coupling region. The usual phonon-induced BCS model is unable 
to produce a large T c in the weak-coupling region. 

Despite much effort, the normal state of the high-T c superconductors has not been satisfacto- 
rily understood. Unlike the conventional superconductors, their normal state exhibits peculiar 
properties. The thermodynamic and electromagnetic observables of these materials above T c 
have temperature dependencies which are very different from those of a Fermi liquid |J. There 
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are controversies about the appropriate microscopic hamiltonian, pairing mechanism, and gap 
parameter for them|| 0. 

The BCS theory considers N electrons of spacing L, interacting via a weak potential of 
short range r such that r « L and r << £. When suitably scaled, most properties of 
the system should be insensitive to the details of the potential and be universal functions of 
the dimensionless variable L/C, ||. In this work we study the weak-coupling BCS problem in 
three dimensions for s, p, and d waves with two objectives in mind. The first is to identify the 
universal nature of the solution appropriate to high-T c superconductors. We would specially 
be interested to find out if the weak-coupling BCS theory can explain some of the universal 
behaviors of high-T c materials independent of the above-mentioned controversies. The second 
objective is to find out to what extent the universal nature of the solution is modified in the 
presence of realistic finite-range (nonlocal separable) potentials. Instead of solving the BCS 
equation on the lattice with appropriate symmetry, we solved the equations in the continuum. 
This procedure should suffice for present objectives. There is also the possibility of Cooper 
pairing in non-s waves, such as, p-wave pairing in superfluid 3 He [0, [H]] and <i-wave pairing in 
some superconductors 0. Hence the present discussion of universality is also extended to p 
and d waves. 

In place of the standard phonon-induced BCS model we employ a renormalized BCS model 
in three dimensions with separable and zero-range potentials, which has certain advantages. 
The standard BCS model yields the following linear correlation between T c and Tp, where Tp, 
is the Debye temperature: T c « 1.1 3Td exp(— 1/A) ||, where A is the effective strength of the 
phonon-induced BCS interaction. Due to the above correlation with Tp>, T c of the standard BCS 
model is low. This correlation between T c and Tp> is fundamental in explaining the observed 
isotope effect in conventional superconductors ||. The high-T c materials exhibit a very reduced 
and negligible isotope effect and a visible linear correlation between T c and Tp ||11||, where Tp 



is the Fermi temperature. The present renormalized BCS model yields a linear scaling between 
T c and T F . Because of this scaling with T F , the present T c can be large and appropriate for the 
high-T c materials in the weak-coupling region. In addition, the present model also produces 
an appropriate T c /Tp ratio and a small £ in the weak-coupling region in accord with recent 
experiments [[1 1] on high-T c materials. 

Previously, there have been studies of the solution of BCS equations in terms of potential 
strength, Vq, or the pair scattering length in vacuum, a, employing a short-range potential [ |T2|| . 
Such studies have not revealed the universal nature of the transition from weak- to medium- 
coupling. Here, we employ the zero-temperature condensation energy per particle, AU, of the 
BCS condensate as the reference variable for studying the problem. As AU increases, one 
passes from weak to medium coupling. We calculate the zero-temperature BCS gap A(0), T c , 
the specific heat per particle C S (T C ) in different partial waves and the zero-temperature pair 
size £ in s wave. These observables obey robust universal scaling as functions of AU valid 
over several decades in the weak-coupling region independent of the range of potential. Similar 
scalings were not found when A(0), T c , C S (T C ) and £ were considered as a function of Vq or a 



as in Ref. 12 



We also calculate the temperature dependencies of different quantities, such as, the BCS gap 
A(T), penetration depth A S (T), spin-susceptibility Xs(T), C S (T), internal energy per particle 
U S (T) and entropy S S (T) for T < T c . Of these, the T dependencies of S S (T), C S (T), Xs(T), 
and A S (T) are interesting. For isotropic s wave, the BCS theory yields exponential dependence 
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on temperature as T — > for these observables independent of space dimension [|, |j~3f . The 



observed power-law dependence on temperature in some of these quantities |L3|, |M[ [T^] can be 
explained with anisotropic gap function in non-s waves with node(s) on the Fermi surface. We 
find universal power-law dependence in non-s waves independent of the range or strength of 
potential. For I ^ we find 

S S {T)/S S (T C ) « (T/T c )^ (1) 
C s (T)/C n (T c ) « D(T/T c f° (2) 
Xs (T)/ Xs (T c ) « (T/T e )* (3) 
AA(T) = (A S (T) - A s (0))/A s (0) ~ (T/T c )^, (4) 

valid for a wide range of temperature. The suffix n and s refer to normal and superconduct- 
ing states, respectively. Similar power-law dependencies were predicted from an analysis of 
experimental data jTJJ] as well as from a calculation based on Eliashberg equation ||15|| . 

From the weak-coupling BCS theory we established the following relations analytically: 
A(0)/ v / A^ = ^8/3, Tjy/AU = ^8/3A-\ G = C S (T C )/^MJ = ^2/3 (tt 2 A' 1 + 1.5AB 2 ), 
£ 2 = A~ 2 (0)/2 = 3/(16A£7), # = (D - 1) = AC/C n (T c ) = 1.5A 2 B 2 /tt 2 , and AU/U n (T c ) 
= 1.5A 2 /tt 2 where AC = C S (T C ) — C n (T c ) and the universal constants A and B are defined 
by A = A(0)/T c and B 2 = -[d{A(T)/A(0)} 2 /d(T/T c )] T=Tc . Unless the units of the variables 
are explicitly mentioned, all energy (momentum) variables are expressed in units of Ep (kp), 
such that /i = fi/Ep, T = T/Tp, q = q/kp, E^ = E^/Ep, E F = kp = k B = 1, etc, where /i is 
the chemical potential and Ep is the Fermi energy. The lengths are expressed in units of kp 1 : 

c = a-,-. 

In Sec. II we derive the present renormalized BCS and number equations. In Sec. Ill we 
present an analytic study of the renormalized BCS equation and find several universal relations 
among the observables. In Sec. IV we present a numerical study of the present equation 
and establish power-law temperature dependence of some of the observables below T c in non-s 
waves. Finally, in Sec. V we present a summary of our findings. 



2 Renormalized BCS and Number Equations 

We consider a weakly attractive short-range potential between electrons in the angular momen- 
tum state Im, 

U pq = — Vogpi m gqi m Yi m {p,p)Yi m {fl q ) , (5) 

where g is the potential form factor and Q (= $0)represents the polar and azimuthal angles. 
This potential leads to Cooper instability for any Vq and Im. In even (odd) partial waves, 
pairing occurs in singlet (triplet) state governed by the Cooper equation 

V, l = E gl m \Y lm {n q )\ 2 {2e q -E)-\ (6) 
q(g>l) 

with B c = (2 — E) the Cooper binding, e q = h 2 q 2 /2m where q is the wave number and m the 
electron mass. 
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At a finite temperature, T, one has the following BCS gap and number equations for iV 
electrons 



A r 



N 



pq 2£ n 2T' 



(7) 



E 



i - 



^ tanh — - § 
2T 



-q - [fa - + 



|Aq| 2 ] 1/2 , 



with A q the 



gap function and /i the chemical potential, 
in a mixed angular momentum state, here 
is formed in a state of well-defined Im, so 



where E, 

Though it is possible to have a BCS condensate 
we assume, as in Ref. [fLOU , that the condensate 
that A q has the following anisotropic form: A q = g q i m A \/4iiYi m (Q q ) where A and g q \ m are 
dimensionless. The BCS gap is defined by A(T) = 
average of A q on the Fermi surface. Equations 
equation 



9q{=i)im^0i which is the root-mean-square 
H) and ([?D lead to the renormalized BCS 
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The summation is evaluated according to 



(10) 



where / dO, = J n d(p Jq sm9d9. Equations @ and @ can be explicitly written as 

16n 
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(12) 



The two terms in Eq. (JTTJ) or (|T^) under integral have ultraviolet divergences. However, the 
difference between these two terms is finite. In the absence of potential form factors (g q i m = 1), 
these equations are completely independent of potential and are governed by the observable 
B c . This is why these equations are called renormalized BCS equations [^, ^6| . The quantity 
B c plays the role of a potential-independent coupling of interaction. 

Now we calculate the critical temperature T c of Eq. (|12|), in the special case g q i m = 1. This 
potential is independent of a range parameter and is usually called a zero-range potential. At 
T = T c , (A(T C ) = 0), Eq. fll2|) can be analytically integrated to yield 



2exp(7 



7T 



2B r » 0.590a 



(13) 
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where 7 = 0.57722... is the Euler constant. If Tp is a few thousand Kelvins, for a small B c 
in the weak-coupling region, one can have T c > 100 K appropriate for some high-T c materials. 
The standard BCS model yields in this case || 

KfZEtoM. (14 ) 

1 D it V Id 

To illustrate the advantage of Eq. ([13]) over (|14]) in predicting a large T c in the weak- 
coupling limit, let us consider a specific example with Tp> = 300 K and Tp = 3000 K. In 
the standard BCS result flUD, the weak-coupling region is usually defined by B c rs 1 meV or 
B c /Td « 0.037. The smallness of -B c justifies the weak-coupling limit and we take S c < 1 meV 
as defining the weak-coupling region. Schreiffer suggested that B c is the proper measure of 
coupling. He noted that B c = 0.1 meV is safely within the weak-coupling domain. In this case 
for B c = 1 meV = 0.0037 one obtains from Eq. ( |I4"D that T c is 46 K. From Eq. (|13"D, we obtain 
T c = 0.036 = 107 K. This result reflects an enhancement of T c in the renormalized model. In 
order to provide further evidence of the weak coupling limit of the present renormalized model 
with T c = 0.036, we solved the number equation ([11]) numerically for the chemical potential /1 
and obtained [i = 1.000, which is in the weak-coupling domain. The renormalized result flT5| ) 
has the advantage of producing the experimentally observed linear scaling between T c and Tp 
in high-T c materials |TT|| . 



3 Analytic Study of the Renormalized BCS Equation 

There is no cut-off in the renormalized BCS equation (|T2"]). At T = Eq. ( |T2"D can be solved 
analytically in the absence of potential form factors: g q \ m = 1. Then each integral in Eq. (|T2| ) 
is divergent at the upper limit A, but for a sufficiently large A the difference becomes finite. 
Now Eq. (|1^) can be integrated in the weak-coupling limit (/i = 1) to yield: 

2y/X - ln(e 2 S c /8) = 2y/X - 2 ln[e 2 A(0) V4?r/8] + In F 2 , 

where InF = - J dtt\Y lm (Q)\ 2 In \Y lm (Q)\ with e = 2.718281... For A -> 00 this leads to A(0) = 
F\l2B c /(ey/ir). However, T c is given by Eq. (p~3|) for all Im. In this case we have the universal 
constant A = A(0)/T c = F0F/{2 exp( 7 )}. 

Though A is independent of interaction model and dimension of space, A(0) and T c are 
dependent on them. For example, for a s-wave zero-range interaction we have A(0) = \/1B c 
and T c = exp('y)y/2B c /7T from a renormalized BCS model in two dimensions0, distinct from 
the above three-dimensional relations. For a fixed coupling, denoted by a B c , we find an 
enhancement of T c in two dimensions over that in three dimensions by a factor of e/2. In both 
two and three dimensions the renormalized BCS equation provides an enhanced T c over the 
standard BCS T c given by Eq. ([14]). 

The entropy of the system is given by 

S(T) = -2£[(1 " U) Ml " /q) + /qln/J, (15) 
q 

where / q = 1/(1 + ex P (EjT)). 
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The condensation energy per particle at T = is given by 



C 2 A 2 

*u = \u.-u n \= £ 2C-£(C 

q(g<i) 



where C 



jj). In the absence of potential form factors this can be evaluated to lead to 



AU=^JdnA 2 (0)\Y lm (m 2 i 



which yields A(0)/v / AZ7 = a/8/3 for all lm. Using the universal relation between A(0) and 



T c , one obtains TjVAU = ^%/3A- 1 . For all lm, U n (T c ) = tt 2 T 2 /A, so that AU/U n (T c 
3A 2 /(2ir 2 ). 

The superconducting specific heat per particle is given by 



C, 



NT 2 



(16) 



The normal specific heat C n is given by Eq. (|1^) with A q = 0. The jump in specific heat per 
particle at T = T c (A(T C ) = 0), AC = [C s - C n ] Tc is given by 



AC 



1 



NT, 
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q l 



/q(l " /q 



dT 



(17) 



In the special case g q i m = 1, the radial integral in Eq. (|17D can be evaluated as in Ref. [@ and 
we get P 



AC 



AT 



dQ 



x 



/q(l " U) 



dA 2 (T) 
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\Y lm (n)\ 2 - 



;is) 



This leads to AC = -(3/4)[dA 2 (T)/dT} T=Tc = (3/A)A 2 B 2 T c for all lm. From this, we obtain 
H = (D - 1) = AC/C n (T c ) = l^A^/vr 2 , where C n (T c ) = tt 2 T c /2, so that C S {T C ) = (vr 2 + 



1.5A 2 B 2 )T c /2. Consequently, C S (T C )/VAU = G = ^J2/3{tt 2 A- 1 + 1.5AB 2 ). The numerical 
values of the constants A, B, H, F and G are given in Table I. 
The spin-susceptibility x °f the system is defined by [T(| 



2/i 



iV 



E/q(Wq), 



(19) 



where /iat is the nuclear magneton. At T = T c , X«CO = Xn(T") and it is appropriate to study 
the ratio x*( T ) I 'Xn{T c ) ■ 

Finally, it is also of interest to study the penetration depth A defined by |§ 



A- 2 (T) = A" 2 (0) 



1 - 



NT 



EA(Wq 



(20) 
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In the numerical study of next section we shall calculate AA(T) = (A(T) — A(0))/A(0). 

The dimensionless s-wave pair radius defined by £ 2 = (ip q \r 2 \ip q ) / {ip q \ip q ), with the pair wave 
function ip q = g q i m A/ (2E q ), can be evaluated by using r 2 = —V 2 . In the weak-coupling limit, 
the zero-range analytic result of Ref. @ leads to £ 2 = A~ 2 (0)/2 = 3/(16A[/). Consequently, 



1 



V2AT C 



(21) 



4 Numerical Study 

Equations ([ll|) and ( |12|) are solved numerically without approximation in s, p and d waves for 
separable potentials with dimensionless form factors g q i m = e l J 2 [a/(e q + a)]^ +2 ^ 2 with correct 



threshold behavior as q — * 0, where a is the range parameter. (Normally, one uses in Eq. ( |10| ) 
de q y/e q = de q ^, ^j.) Following Refs. 0, [UJ, we calculated A(0), T c , C,s(T c ), the s-wave pair 
radius £ 2 at T = as well as A(T), A(T), C(T), S(T), and U(T) for different Vo and a. In Fig. 
1 we plot A(0), T c , C S (T C ), and £ 2 versus At/ and establish universal scalings mentioned before. 
The calculations were repeated by varying a from 1 to oo and we found Fig. 1 to be insensitive 
to this variation for each Im. For / ^ 0, Eqs. flTT| ) and ( |T2"D diverge for a — > oo and calculations 
were performed for a = 1 to 10. The increase in AU of Fig. 1 corresponds to an increase in 
coupling V . We could plot the variables of Fig. 1 in terms of V as in Ref. Then each 

a leads to a distinct curve. However, if we express the variation in Vo by a variation of an 
observable of the superconductor, such as AU or T c , universal potential-independent scalings 
are obtained. In each case the exponent and prefactor of each scaling relation are in excellent 
agreement with the analytic result obtained above without form factors. 

TABLE I. Numerical values of various constants and exponents in different angular momen- 
tum states 
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1.764 
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1.43 


11.11 
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2.8563 
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2 
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3.3300 
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10.59 
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2.6 




2.4 


20 


2.7748 


1.382 


1.57 


0.72 


10.00 
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1.2 


1.1 


21 


3.1006 


1.544 


1.63 


0.96 


10.24 


2.1 


2 


1.4 


1.5 


22 


3.1006 


1.544 


1.63 


0.96 


10.24 


2.1 


2 


1.4 


1.5 



The T c should not arbitrarily increase with coupling as Fig. 1 may imply. With increased 
coupling the electron pairs form composite bosons which undergo Bose condensation below 
T = T c = 0.218, for bosons with twice the electron mass and half the electron density 
Hence the T c curve of Fig. 1 is only plotted up to about T c = 0.1. For a large class of 
unconventional three-dimensional superconductors T c has been estimated to be 0.05 [11], where 
the universality of the present study should hold. For a typical high-T c material T c = 0.04 and 
from Eq. ([H]) we find pair-size £ ~ 10 in s wave. Hence with the increase of T c , £ has reduced 
appropriately in the weak coupling region as found experimentally. 

Next we studied the temperature dependence of A(T), S S (T), C S (T), Xs(T), U S (T), and 
A S (T) for T < T c for different Vo, and range a varying from 1 to oo. We found that A(T)/A(0) 
versus T/T c is an universal function for each Im independent of potential parameters. We 
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find the universal fit A(T)/A(0) ~ B(l — T /T c ) 1 ^ 2 valid for T pa T c with numerical values of 
B quoted in Table I. For s-wave BCS superconductors S S (T), C S (T), Xs(T), and A S (T) have 
exponential dependence on T as T — > 0. But for non-s wave states, these variables have power- 
law dependence on T as observed in some materials |13|, [14], [15]. In Figs. 2, 3, 4, and 5 we plot 
S s (T)/S n (T c ), C s (T)/C n (T c ), Xs(T)/xn(T c ), and AA(T)/A(0), respectively, versus T/T c where 
AA(T) = (A(T) — A(0))/A(0). As commented in Ref. []nj, \s will be significantly different 
from Xn only for even I. We have calculated Xs only for I = and 2. Scalings ([!]) — (|j) are 
established in Figs. 2 — 5. The exponents of these scalings are given in Table I. In order to 
find (3c we also plotted C s (T)/C n (T c ) versus T/T c on log scale. That plot was used to calculate 
the exponent (3c- However, on log scale different curves nearly overlap and hence that plot is 
not shown here. From Fig. 3 we find that the zero of [C S (T) — C n (T)] appears in all cases for 
T/T c pa 0.5. Moreover, all curves for superconducting specific heat meet at T/T c pa 0.6. These 
two behaviors seems to be typical for models based on BCS equations. 

The constants in Table I for different Ira satisfy C 00 > C\\ > C 2 \ = C22 > C10 > C20, where 
C\ m stands for F, A, B, H, and G. Hence the following sequence of Im states represents the 
increase of anisotropy for the gap function: 00, 11, (21,22), 10, and 20. From the plot of 
entropy in Fig. 2, we find that this sequence of Im also represents the increase of disorder 
and consequently, a decrease in superconductivity or an approximation to the normal state, 
as is clear from Figs. 3 — 5. Because of approximation to more anisotropy and disorder, the 
observables for the normal state are closer to the superconducting I 7^ state than to the 
superconducting I = state. 



The exponents (3s, (3c-, P\ anc ^ /^A are critical exponents near T = T c . Wilson |L8| discussed 
the universal nature of similar critical exponents in ferromagnetism and concluded that the 
numerical value of those exponents depend on the dimensionality of space and the symmetry of 
the order parameter of phase transition. Recently, we have calculated some of these exponents 
using the renormalized BCS equation in two dimensions [I7[] . From these studies it seems 



that these universal critical exponents of BCS superconductivity are also determined by the 
dimensionality of space and the symmetry of the order parameter A q . 



5 Conclusion 

Through a numerical study of the renormalized weak-coupling BCS equation in three dimen- 
sions in s, p and d waves we have established robust scaling of A(0), T c , C S (T C ), and £ 2 , as 
a function of AU, independent of range of a general separable potential. The T dependence 
of S S (T), C S (T), Xs{T), and AA(T) below T c in non-s waves show power-law scalings distinct 
to some high-T c materials 0, 0, |15]| . No power-law T dependence is found in s wave for 
these observables. The universal nature of the solution does not essentially change with the 
potential range and remains valid for a zero-range potential. In the weak-coupling limit the 
present solutions of the renormalized BCS equations simulates typical high-T c values for the 
coherence length £, and T c . They also exhibit the T c versus Tp linear correlation (at a fixed 
B c ) as observed by Uemura [|ll[ . 

We thank John Simon Guggenheim Memorial Foundation, Conselho Nacional de Desen- 
volvimento Cientifico e Tecnologico and Fundagao de Amparo a Pesquisa do Estado de Sao 
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Figure Captions: 

1. C S {T C ) (dashed line), T c (dotted line), A(0) (dashed-dotted line) for different Im and 
s-wave pair radius £ 2 (solid line) versus zero-temperature condensation energy A?7 for different 
Vq and a (from 1 to oo). For C S (T C ) and T c there are six distinct lines and for A(0) we have 
a single line for all a and Im. The lines for C S (T C ) (T c ) correspond to Im = 00, 11, (21, 22), 10, 
and 20 from top to bottom (bottom to top). 

2. Entropy S S (T)/S S (T C ) versus T/T c for different Im, Vq, and a between 1 to oo. The 
curves are labelled by Im. 

3. Same as Fig. 2 for specific heat C s {T)/C n {T c ) versus T/T c . 

4. Same as Fig. 2 for spin-susceptibility Xs(T)/Xs(T c ) versus T/T c . 

5. Same as Fig. 2 for AA(T) versus T/T c . 
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